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We set up a melting model for vortex lattices in high-temperature superconductors based on 
the continuum elasticity theory of the vortex lattice. The model is Gaussian and includes defect 
fluctuations by means of a discrete-valued vortex gauge field. We derive the melting temperature of 
the lattice and predict the size of the Lindemann number. Our result agrees well with experiments 
on YBa2 Cus Or-i, but not on Bi2 Sr2 CaCu208. We also calculate the jumps in the entropy and 
the magnetic induction at the melting transition. 
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rson? ^^^^ lattice melting, proposed by Nelson in 1988 
Uj , has become an important topic in the phenomenology 
of high-temperature superconductors. The properties of 
this transition have been studied in a large number of 
theoretical as well as experimentel papers. The simplest 
method to estimate the temperature where the transi- 
tion takes place is adapted from the famous Lin^eiftfinn 
criterion of three dimensiong.|^o]|'dinary crystals IBJ . The 
adaption by Houghton [2] to vortex lattices states that 
the vortex lattice undergoes a melting transition once 
the mean thermal displacement (u^)^/^ becomes a certain 
fraction of the lattice spacing a « (<I>o/i3)^/^ where $0 is 
the flux quantum, B is the induction field. The size of the 
Lindemann number cl — (u^)^/^/a is not predicted by 
Lindemann's criterion and must be extracted from exper- 
iments, and is usuajjvfo^^ in the range ~ 0.1—0.3. It 
has been shown in 12] that cl is, in fact, rougly universal, 
with only a weak dependence on the induction of the sys- 
tem. The most prominent examples of high-temperature 
superconductors exhibiting vortex lattice melting are the 
anisotropic compound YBa2 Cua 07-5 (YBCO), and the 
strongly layered compound Bi2 Sr2 C^,Cji208 (BSCCO). 
Decoration experiments on BSCCO f'sfshow the forma- 
tion of a triangular vortexlatti^e, on YBCO of a tilted 
square lattice of vortices 101 close to the melting region, 
the latter being favored by the c?-wave symmetry of the 
order parameter l'5[7* An explicit calculation of the Lin- 
demimji|.mj|mber = {"^^^l l^-l^^or YBCO can be found 
in (liflsfTor BSCCO in f^r^ 

In this letter, we calculate the size of the Lindemann 
number cl from a Gaussian model which takes into ac- 
count the lattice elasticity and the defect degrees of free- 
dom in the simplest possible way. The relevance of defect 
fluctiiaj^kj.^Sj^to vortex melting has also been emphasized 
by i r^J. lor ordinsu^-gj^ystals, this is meanwhile textbook 
material: in Ref. E one of us has set up simple Gaussian 
lattice models based on linear elasticity and a fluctuat- 
ing discrete-valued defect gauge field for square lattice 
crystals which clearly display melting transitions. These 
are usually of first order but, if rotational stiffness is suffi- 
ciently high, they may be split into two successive contin- 
uous transitions. An important virtue of these models is 
that in the first-order case, where fluctuations are small. 
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FIG. 1: Melting curve B = B„(r) fofe^ ^BgQ^Mi d BSCCO. 
iTJ:^^^xperimental values are from Ref. Il7 tor VBCO and Ref. 
ibO tor BSCCO. The numbers at the vertices of the making 
curve are the Lindemann numbers cl calculated from 

they lead to a simple universal melting formula determing 
the melting point in terms of the elastic constants. The 
universal result is found from a lowest-order approxima- 
tion, in which one identifies the melting point with the 
intersection of the high-temperature expansion of the free 
energy density dominated by defect fluctuations with the 
low-temperature expansion dominated by elastic fluctu- 
ations. The resulting universal formula for the melting 
temperature determines also the size of the Lindemann 
number. In two dimensions, this procedure was„rc(^cpjly 
generalized from square to triangular lattices ISJ. ihe 
intersection criterium was also used before to find the 
melting point of vortex lattices using the AbrikpsoVj^ ap- 
proximation of the Ginzburg-Landau model Llu] useful 
for YBCO. Here we shall apply our defect model to cal- 
culate the melting curve of the vortex lattices in YBCO 
and BSCCO. 

Due to the large penetration depth Xab in the layers 
in comparison to a we have to take into account the full 
non-local elasticity constants when integrating QXP^^dily 
Fourier space, as emphasized by Brandt in Ref. llT. i'or 
our Gaussian model, the partition function of the vor- 
tex lattice can be split into Z — ZqZ-a where Zq is the 
partition function of the rigid lattice and Zi^ is thermally 
fluctuating part calculated via the elastic Hamiltonian 
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plus defects. Due to the translational invarince ofhtiiedie Inserting the defect gauge field, it is extended to * 

vortex lattice in the direction of the vortices, whichtw#i 

shall take to be thi^^^^^axis, we may simply extend"^?!?^" 

moddg j^n^jquare \%\ and triangular and triangular ^Ig^- 

tices fjSj by ^ third dimension along the z-axis, which we 

artificially discretize to have a lattice spacing a^, whose 

value will be fixed later. The elastic energy is 

-Bel = - ^ J dz [(cii - 2c66)(ViWi)^ + ceeC^iUj + VjUi)/2 + c^^id^Ui)'^] . 
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In the canonical stress representation, the partition func- tion becomes 
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The subscripts i, j (/, m, n) have the values 1, 2 (1, . . . , 3), 
the the-dimcnsional vactors Ui (x) are give the transverse 
displacements of the line elements of the vortex lines 
with coordinate x, and crim(x) is the conjugate stress 
fields. We have suppressed the spatial arguments of the 
elasticity parameters, which are really functional matri- 
ces Cy (x, x') = Cy (x — x'). In momentum space, the 
brackets are diagonaL-^jj^^cir precise forms were cal- 
culated by Brandt ITTj. The functional matrix /3 is 
given hy f3 = v ce6/kBT{2Tr)'^ where v is the volume 
0^03 (square) or a^a3^/3/2 (triangular) of the funda- 
mentaljXell. For a square lattice, the lattice derivates 
Vi in ^ are given by Vi/(x) = [/(x -I- ae^) — /(x)]/a 
and V3/(x) = [/(x + 0363) - /(x)]/a3. For a triangu- 
lar lattice, the xy-part of the lattice has the link vec- 
tors ae(m) with 6(13) = (cos 27r/6, ± sin 27r/6, 0) and 
6(2) = (~ljO, 0). The lattice derivates around a plaque- 
tte are defined by V(i)/(x) = [/(x + ae(i)) - /(x)] /a. 



V(2)/(x) - [/(x)-/(x-ae(2))]/a, V(3)/(x) - 
[/(x — ae(2)) — /(x -|- ae(i))] /c. From these we define 
discretB, cartesian derivates used in the partition func- 
tion m V./(x) = (2/3)e(,),V(i)/(x) and V3/(x) = 
[/(x + a3e3) - /(x)]/a3.* check 

We, are now prepared to specify the matrix Dij (x) in 
Eq. It is a discrete-valued local defect matrix com- 
posed of integer- valued defect gauge fields ni,n2,n2 for 
square and triangular vortex lattices as follows: 

^^^^ 

The vortex gauge fields sufaces in 

units of the Burgers vectors 10, Vf. By summing over 
all rii,2,3(x), the partition function includes all defect 
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fluctuations, dislocations as well as disclinations. There 
is a constraint for a vortex lattice which does not exist 
for ordinary three-dimensional lattices. Dislocations in 
the vortex lattice can be both screw or edge type, but 
in cither case the defect lines arc confined in the plane 
spjiiingd ]!jg^(Jjjj^!fiPurger's vector and the magnetic field 
[|iT2riT3i|.' The reason is that the flux lines in a vortex 
lattice cannot be broken. This results in the constraint 
Dii = —D22 on the defect flelds. 

We now calculate the low-temperature and high- 
temperature expansions of the partition function Zf\ to 
lowest order, which includes only the — 0-term. By 
carrying cnit. the integration cxy^,4:l^^(^ijjl|acement fields 
Mi(x) in Igf we obtain, as in the leading term in 

the low-temperature expansion of the free energy 
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HereJfm is the eigenvalue of iVm- The k, fcs-integrations 
in run over the BriouUin zone of the vortex lattice of 
volume Vbz = (27r)^/u, as indicated by the subscript BZ. 

Next we calculate the high-temperature expansion 
^T^oc lowest order. By carrying out the integration 
over the displacement fields Wi(x) in ^ and further by 
summing over the defect fields under the considera- 
tion Dii = —D22 mentioned above we obtain 
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The constat C has the values Co = 1 for the square 
vortex lattice and Ca = VS for, the triangular lattice. 

J atrain ° 

From the partition function (0) with no defects (rim = 



0) we obtain for the Lindemann number cj^ 
the momentum integral 
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This expression can be simplified by takin&,_int^2^c- 
count that cn is much larger than C66,C44 
shall calculate the melting temperature from the intersec- 
tion of low-and high-temperature expansions, obtained 
by equating Z'^'^^ — Zl^°° . By taking into account 
det[a|V3V3] = 1 we obtain,.^, In <C lea a-nd further that 
the summand for i = 1 in fiuj) is much smaller than the 
summand for i = 6. In the following analytic discussion 



(but not in the numerical plots) we will neglect the cor- 
responding terms. The temperature of melting is then 
given by 
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where det[c66] is the determinant of the N x N functional 
matrix cge- The elastic moduli C44 and cge at low reduced 
magnetic fields Jt^^^^ / Hc2 < 0.25 can be taken from 
Brandt's paper iTff 
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where Ac is the penetration depth in the a;y-planc, and 
C = 1. Kbz is the boundary of the circular Brillouin 
— zone K'^,^ = ATTB/(f)o- At high fields (5 > 0.5), cge is 
.i^altered by the factor C k. 0.71(1 — 6), and the penetration 
depths in cee, C44 are replaced by = X"^ /{\ — b), where 
A denotes either Xab or Ac. In addition, the last two 
termsLQf Cji^ are replaced by B(t)Q/lQi{^X^^. For YBCO SfCjj^jjaj^]^ 
have mf^{T) = A(0)[1 - (^/Tc)]-l/^ for BSCCO I'^f" ^ 
A(0)[1 - (T/rc)2]"i. For the calculation of C44 in 
we have used a momentum cutoff in the two-vortex 
interaction potential k < 2/ (u^)^/"^ * andiix>t^|y inverse 
of the correlation length l/£,^m Ref. IHir I'he cutoff 
is due to thermal softening IB], and becomes relevant for 
SO^'^y/'^/C > 1, or equivalently for clV2jt y/H^^^jT)jB > 
1, which is fulfilled in the melting regime of BSCCO, but 
not for YBCO. 

It remains to determine the effective lattice spacmg 
03 along the vortex lines. The melting condition 1)111 |l 
erges for 03 — > 0, whereas the displacement average 
remains finite in this limit. An elementary defect 
arises from an interchange of two vortex strings in the 
sample. This takes place over a typical length scale in 
the z-direction. To properly count the defect degrees of 
freedom in the partition function, the length 03 should 
coincide with this length scale. The length ap is deter- 
mined from the condition that the sum of elastic dis- 
placement energy and the energy required ^"JJijf ^tch the 
line against the line tension is minimal ihus we in- 
sert the variational ansatz for the transverse displace- 
I 70i| ient Ui = (5i^i^o exp[— 2|z|/a3] into the elastic en- 
ergy and approximate — V2 ~ (^2) ~ -^o/4 in 
and « (^^) ~ -^o/^ in the elastic constants, where 
the average (•) was taken with respect to a circular Bril- 
louin zone. The optimal length Oz is chosen such that 
Ec\ is minimal for a fixed amplitude ~ a.* 

In the following, we treat first the more isotropic 
square vortex crystal YBCO (a — y^^o/S). From 
and C44 for YBCO, the optimal length is given by 
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CI3 = 4aoAab/AcC(l — ^)^^^-* When coragaring the melting 
criterium of the defect model in Eq. with the Linde- 
mann criterium obtained by equating the parameter ^ro\i 
to a fixed number, we get almost a precise coincidence of 
the two criteria. An exact coincidence is mached when 
taking into account that the integrand in flU)) and in /gg 
of Eq. (|7|' receive their main contribution from the re- 
gion k w {k?) w i^Bz/v^ where we can put k. 0, 
resulting j|t^ 03056/^^044 ~ 2. Using this approximation 
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whd;ihe length of the dislocation. Remebering this we ob- 
is taila by a similar procedure as for YBCO the disloca- 
length 03 w Aa{^/2iT)\ab / ^c- From this we find 
oIcqq/o^caa ^ 1, resulting in Zge ~ 0. By taking into 
account that i?7r^A^/32(/)o <C 1 we obtain that C/^Ajk Jfyi ) 
for \k^\ < tt/c is dominated by the last term in 
Then we obtain * 



in Eqs. fj' and I|I1UI) . we can perform the integrals nxi- 
mcrically, adnm obtain from the melting criterium 
preciselyite Lindemann criterium with the Lindcmann 
number fiu)) fixed to be 
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Denoting by the spacing between the Cu02 layers we 
obtain for the entropy jump per layer and vortex 



cl 



ksTyn ■ 0.036 



fcsT^c- 1.60 



C66(^,0)C44(^,0) 



Km 



ASd « kBTrn{d/dTrn)ias/as) \n[Z^ 
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Finally, we make use of the Clausius-Clapeyron equation 
ASd/a^as = -{dHm/dT)AB/A'K, relating the jump of 
the entropy density ASd/{a?s) to the jump of the mag- Bm{T) 
netic induction AB across the melting transition. Here 
H„i is the external magnetic field on the melting lir 
FronX|thc the Clausius-Clapeyron equation and Eqs. 
and fibjl . we obtain the following relations near Tc (where 
b is small*) 



With the help of our melting criterion nlB|l and 
the Clausius-Clapeyron equation {weMse^dH,n/dT w 
m^n/dT based on B^{T) > H^iiT) ^^^fvte obtain for 
BSCCO; 
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parameter cl, as well as the entropy ,and mag 



These results agree with the scaling results in Ref. l(TSj7" 
where multiplicative prefactors are determined by fits to 
experimental curves (there is only a slight discrepancy 
because we use a different temperature dependence of 
the penetration depth). 

Next, we calculate the corresponding expressions in 
the case of the more layered crystal BSCCO (a — 
3i/4^2i/2y/0j^/iJ). First, we have to determine the dis- 
location length 03 in this case. For dislocation moves 
we have (u^)^/^ ^ l/K^z* This n^ns that we can ne- 
glect the last two terms of C44 in Ijllal) . coming from the 
self-energy of the vortex line, for the determination of 
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loB/jUging these values we obtain by numerical integration 
of (fID|) 
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The first term comes from the integration region jfesj < 
jf), the second from the region 1/Xab < l^sl < tt/c in 
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From the experimental values Aafc(O) ~ 1500A 
(Aq6(0) « 22OOA), the Cu02 layer spacing as = 12k 
(s = I4A), Tc = 92.7K (Tc = 90K) and the anisotropy 
parameter 7 = Xc/Xab ~ 5 (j — Xc/Xab ~ 200) for 
optimally doped YBCO (BSCCO) we now calculate ex- 
Ijicit the melting curves with the associated Lindemann 

tic in- 
. The 



duction jumps ASd and AB from hUyii and 
latter are found from the intersection criteriumj^ the 
the full low- and high-temperature expressions ^ and 
(^,respectively. The results are shown in Figs, l^pd 2, 
rgspectively. Our approximate analytic results fiV^ and 
fitjl) turn out to give practically the same curves. For 
comparisons, we show vq, ^IjAtJr -fi g^yjj^jf ^jfiig^g^j^ ^ntal 
£Ufy.es^f(^^^]g.Cp of Ref . Ilf ; TS ^ tfj and^br BSCCO of Ref. 
120721. ' I he temperature dependence of the Lindemann 
parameter for BSCCO seen in Fig.l comes mainly from 
the second term in . which is the largest near T^. Our 
curves are in reasonable agreement with with the exper- 
imental data for YBCO but not for BSCCO. The smaU 
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FIG. 2: Entropy jump per layer per vortex ASd (first row) 
and jump of magnetic induction field AS (second row) at the 
melting tr^j^^^ij^j^ The experirrk(;j[i|tej values for YBCO are 
from Ref. jiV tOiCjr^^^j^Jjd R-ef- 1X8 for AS by squid experi- 
ments and Ref. liy by torque measurop^n^^The experimi 
tal values for BSCCO are from Ref. LjuTl squid) and Ref. 
(torque). 

discrepancies for YBCO can be explained by the simplic- 



ity of our model and the roughness of the approximations. 
The largest discrepancy of the melting curves for BSCCO 
lies in the high-field regime. Due to the large anisotropy 
of BSCCO, the layers decouple in this low-temperature 
rgrime^nto two dimensional l^tiGGS 9^ (Jl%5manl^^ vortices 
l|!^3f bv Josephson decouphng iBrf^r "J-l^i^ process is ig- 
nored in the simple vortex lattice model set up in this 
paper. In the case of the entropy jump ASd and mag- 
netic field jump AB we obtain in the low-temperature 
regime a rough agreemejjt^of ^oi^r theory with the curves 
of Kadowki .jet.al. ,Rcf. |2T but not with those of Zeldov 
et al. Ref. tTD. The discrepancy iPfjthc^Jwo experimen- 
tal curves may be due to pinning l,T5f. i^br temperatures 
near Tc, the experiments give large values of ASd and of 
AB, which are not reproduced by our model where ASd 
is almost constant at large temperatures. This must be 
due to system degrees of freedom not included in our vor- 
tex lattice model coming from thermally created vortex 
loops in addition to the magnetically created ones form- 
ing the lattice. These vortex loops screen the effective 
repulsive interaction between the magnetic induced .yor- 
owkiF^^ resulting in a rise of the magnetic jump AB 
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We set up a melting model for vortex lattices in high-temperature superconductors based on the 
continuum elasticity theory. The model is Gaussian and includes defect fluctuations by means of 
a discrete-valued vortex gauge fleld. We derive the melting temperature of the lattice and predict 
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soff ^^"^^ lattice melting, proposed by Nelson in 1988 
has become an important topic in the phenomenol- 
ogy of high-temperature superconductors. The proper- 
ties of this transition have been studied in a large num- 
ber of theoretical as well as experimentel papers. The 
simplest method to estimate the temperature where the 
transition takes place is adapted from the famous Lin- 
jdg^^nn criterion of three dimensional ordinary crygta^^^^^ 
IQ. The adaption by Houghton et al. and Brandt Vi\ to 
vortex lattices states that the vortex lattice undergoes a 
melting transition once the mean thermal displacement 
(m^)^/^ becomes a certain fraction of the lattice spac- 
ing a K, ($0/^)^^^ where $0 is the flux quantum, B 
is the induction field. The size of the Lindemann num- 
ber Ch — {v?')^/'^ /a which should be independent of the 
magnetic field is not predicted by Lindemann's criterion. 
It must be extracted from experiments, and is usually 
found in the range w 0.1 — 0.3. The most promi- 
nent examples of high-temperature superconductors ex- 
hibiting vortex lattice melting are the anisotropic com- 
pound YBa2Cu307-_5 (YBCO), and the strongly layered 
compound Bi2Sr2Ca/J^u^08 (BSCCO). Decoration exper- 
iments on BSCCO I'sf^how the formation of a triangu- 
lar vortex^lattjpe, on YBCO of a tilted square lattice of 
vortices llj close to the melting region, the latter being 
favored by the d-wave symmetry oLt^e order parameter 
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FIG. 1: Melting curve B = Bm{T) fo^^ 



iXti^^xperimental values are from Ref. 
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and the anisotropy of the crystal iBj!"* An explicit cal- 
culation of the Lindemann Byif^t*bnf^ ^ (u^)^/^/a fft^ 
YBCO can be found in Ref. |,''a5dfor BSCCO in Ref. 4" 



In this letter, we calculate the size of the Lindemann 
number cl from a Gaussian model which takes into ac- 
count the lattice elasticity and the defect degrees of free- 
dom in the simplest possible way. The relevance of defect 
fluctiiaj^kj.^Sj^to vortex melting has also been emphasized 
by lT2j. For ordin^j^jg'ystals, this is meanwhile textbook 
material: in Ref. E one of us has set up simple Gaussian 
lattice models based on linear elasticity and a fluctuating 
discrete- valued defect gauge field for square lattice crys- 
tals which clearly display first order melting transitions in 
three dimensions. An important virtue of these models is 
that in the first-order case, where fluctuations are small, 
they lead to a simple universal melting formula determing 



are the Lindemann numbers cl calculated from lEf. 

the melting point in terms of the elastic constants. The 
universal result is found from a lowest-order approxima- 
tion, in which one identifies the melting point with the 
intersection of the high-temperature expansion of the free 
energy density dominated by defect fluctuations with the 
low-temperature expansion dominated by elastic fluctu- 
ations. The resulting universal formula for the melting 
temperature determines also the size of the Lindemann 
number. In two dimensions, this procedure was„rc(^ej|;ly 
generalized from square to triangular lattices vd\. ihe 
"Trrtersection criterium was also used before to find the 
melting point of vortex lattices using the AbrikpsoVj^ ap- 
proximation of the Ginzburg-Landau model fi'luif^useful 
for YBCO. Here we shall apply our defect model to cal- 
culate the melting curve of the vortex lattices in YBCO 
and BSCCO. 

Due to the large penetration depth Aafc in the layers 
in comparison to a we have to take into account the full 
non-local elasticity constants when integrating QKyjjjily 
Fourier space, as emphasized by Brandt in Ref. Ifl. i'or 
our Gaussian model, the partition function of the vortex 
lattice can be split into Z — Z^Zf^ where Zq is the parti- 
tion function of the rigid lattice and Zf^ is thermally fluc- 
tuating part calculated via the elastic Hamiltonian plus 
defects. Due to the translational invariance of the vor- 



2 



tex lattice in the direction of the vortices, which we shall 
take to be.g^g2-^-axis, we may simplj(^ e^t^jpd the models 
on square Ip] and triangular lattices raj by a third dimen- 
sion along the z-axis, which we artificially discretize to 
have a lattice spacing 03, whose value will be fixed later. 
The elastic energy is 

■E^ci = ^ ^(ViuO(cii - 2c66)(ViUi) (1) 

X 

+ + VjmO C66 (ViUj + VjU.i) + {dzUi) C44 {dzUi). 

Th subscripts «,j {l,m,n) have values 1,2 (1,...,3), 
the vectors Ui(x) are given by the transverse displace- 
ments of the line elements of the vortex lines with co- 
ordinate X. We have suppressed the spatial arguments 
of the elasticity parameters, which are really functional 
matrices Cy (x, x') = Cij(x p_xQ^2 Their precise forms 
were calculated by Brandt IttJt v is the volume a^aa 
(square) or 0^03-^/3/2 (triangular) of the fundamental 

Sll. For a square lattice, the lattice derivates in 
I are given by Vi/(x) = [/(x + ae,;) — /(x)]/a and 



V3/(x) = [/(x -I- 0363) - /(x)]/a3. For a triangu- 
lar lattice, the a;i/-part of the lattice has the link vec- 
tors ae(„i) with 6(13) = (cos 27r/6, ± sin27r/6, 0) and 
6(2) — (— IjO, 0). The lattice derivates around a plaque- 
tte are defined by V(i)/(x) = [/(x + ae(i)) - /(x)] /a, 
V(2)/(x) = [/(x)-/(x-ae(2))]/a, V(3)/(x) = 
,[/(x — ae(2)) — /(x + ae(i))] /a. From these we de- 
Ifine discrete cartesian derivates used in the Hamilto- 
nian O V,/(x) - (2/3)e(,),V(,)/(x) and V3/(x) = 
[/(x + a3e3) — /(x)]/a3 transforming like the continuum 
deri'ratc^s ^\yith respect to the symmetry group of the lat- 
tice [9]. Therefore, the Hamiltonian has the full sym- 
metry of the triangular lattice and the correct continuum 
elastic energy for zero lattice spacing. 

Within the elastic approximation the displacement 
fields are restricted to values within the fundamental cell. 
In order to contain also defect degrees of freedoii^^^j^ige 
has to put in ^ integer valued defect gauge fields l§J. in 
the canonical stress representation, the partition function 
containing these fields becomes 



Zfi = det 



C66 



4(cii - Cm) 



1/2 
det 



27r/3 



5/2 
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i<m 



n E 

m n„,(x) = -c 





" du' 








exp 1 
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The functional matrix /3 is.^jfj^ by /3 = w cm/kBT{2'n:Y . 
aij represent stress fields \^ . We are now prepared to 
specify the matrix Dij{x) in Eq. (j^. It is a discrete- 
valued local defect matrix composed of.jflJj;j|;cr- valued de- 
fect gauge fielc ^. ^t^,^2 , for square I^Jahd triangular 
vortex lattices as follows: 



ni 
"3 



1 

V3 



n2 



(711-712) + ^rt3 



1 

V3 



{■ni~n2) + -^n3 
-~\rii — 7i2 



(4) 



We now calculate the low-temperature expansion of 
the partition function Zf[ to lowest order, which includes 
only the rim = 0-term. By carrying out the integration 
py^^t|jc^^i||)lacement fields Ui (x) in ^ we obtain, as in 
IS""2.L' the leading term in the low-temperature expansion 
of the free energy 



(3) ^ 



2N 



1 



20 



a J det [(27r/3)c44/c66] 



-NY. 



£{1,6} 



la 



(5) 



The vortex gauge fields specify the Volterra sufaces in 
units of the Burgers vectors. By summing over all 
'^i,2,3(x), the partition function Zn includes all defect 
fiuctuations, dislocations as well as disclinations. There 
is a constraint for a vortex lattice which does not exist 
for ordinary three-dimensional lattices. Dislocations in 
the vortex lattice can be both screw or edge type, but 
in either case the defect lines are confined in the plane 

Pui'^S^^ fecfeWiP^'^S^"'^'^ vector and the magnetic field 
rT3|.' I'hc reason is that the flux lines in a vortex 
lattice cannot be broken. This results in the constraint 
£'11 = —D22 on the defect fields. 



where 



1 1 /■ n 

^" = nTT- / « fc*3 log 
^ VBZ Jbz 



C44 



(6) 

HergJfm is the eigenvalue of iVm- The k, fc3-integrations 
in run over the BriouUin zone of the vortex lattice of 
volume Vbz — {2n)^/v, as indicated by the subscript BZ. 

Next we calculate the high-temperature expansion 
^T^oo lowest order. By carrying out the integration 
over the displacement fields Ui(x) in ^ and further by 
summing over the defect fields 71^ under the considera- 
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with 



-Di2 mentioned above we obtain 



Ka 



N 



2^ det [(2^/3)2c44/c66] 



1 1 /■ 2 

h = --— / d kdkz log 
^ Vbz Jbz 



1 + 



Cll - C66 K*Kj 



C44 KSK3 



(7) 



(8) 



The constant C has the values = 1 for the square 
vortex lattice and Ca = forAhe triangular one. 

From the partition function ^ with no defects {um = 
0) we obtain for the Lindcmann number cj^ = {u'^)^^'^ /a 
the momentum integral 



a?v Vbz 



2 1 
d kdk^ 

BZ C44 



1 



i=l,6 ■ 



C44 J ■> 



(9) 
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This expression can be simplified by takin&,_ir 
count that Cn is much larger than C66,C44 ITI 
shall calculate the melting temperature from the intersec- 
tion of low-and high-temperature expansions, obtained 
by equating Z^^'^ = By taking into account 

det[a|V3V3] = 1 we obtain igZn ^ Zgg and further that 
the summand for z = 1 in Q is much smaller than the 
summand for i = 6. In the following analytic discussion 
(but not in the numerical plots) we will neglect the cor- 
responding terms. The temperature of melting is then 
given by 



knT 



1 



det^/^ [c66 



C = 



(10) 



where det [cee] is the determinant of the N x N functional 
matrix cgg. The elastic moduli C44 and cge at low reduced 
magnetic fields ^^= ^^S /Hc2 < 0.25 can be taken from 
Brandt's paper 1(1 ij 



C66 



C44 



60C 



(STrAah) 



(11) 



In- 



2Xi 



'An{l + Xlk-^ + \l^kl) 32^2^2 i + xiK] 



BZ 



32n^\t,ki 



In 



1 + Ag,fc32/(1 + A2,j^2^) 



(12) 



where Ac is the penetration depth in the xy-plane, and 
C = 1. Kbz is the boundary of the circular Brillouin zone 
K^.^ — 4:TrB/(j)Q. At high fields {b > 0.5), cge is altered 
by the factor C w 0.71(1 — b), and the penetration depths 
in C66, C44 are replaced by A2 = A2/(1 — 6), where A de- 
notes either Aafc or Ac. In addition, the last two terms 
f C4A are replaced by i3(/)o/167r2A2. For YBCO we have 

^xn = A(o)[i - iWLi^'^ ^WT)) = mil - 



(T/rc)]-i/2 forBSCCOl|pf(T7== A(0)[l-(T/rc)4]-i/2^ 

C(r) = i|(p)[l - (r/rc)4]-i/2 . For the calculation of 
C44 in hIt2|I we have used a momentum cutoff in the 
two-vortex interaction potential k < 2/{u'^y^^ * and 



Bpt^^^ inverse of the correlation length l/j^irgg2''^ 
In. ihe cutoff is due to thermal softening 
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riomes relevant for 



Ref. 

and be- 
^ Ij or equivalently for 
■ 1, which is fulfilled in the melt- 



ing regime of BSCCO, but not for YBCO. 

It remains to determine the effective lattice spacing 03 

along the vortex lines. An elementary defect in the 

^^irortex lattice (for example an interchanging of two vor- 
tex strings) takes place over a typical length scale in the 
^-direction which is determined from the condition that 
the sum of elastic displacement energy and the energy 
required to stretch the line against the line tension is 
minimal. A dislocation network containing this length 
spacing in z-direction has to take into account in the 
continuum elastic Hamiltonian for small lattice displace- 
ments. This dislocation network consists of integer fields 
coupling to the displacement fields of length scale larger 
:^r equal to the length spacing of the dislocation network, 
""^he displacement fields of smaller length scales can be 
easily integrated out giving an additive correction term 
to the free energy not relevant for the melting transi- 
tion properties. The relevant part of the free energy 
given by the discretisized free energy — log{ZoZfi) / f3 
with 03 given by the length spacing of the dislocation 
network in z direction. To determine a^, we insert the 
variational ansatz for the transverse displacement field 
Ui — Ji^i^o exp[— 2|z|/a3] into the ccgitinuum version (in 
z direction) of the elastic energy Igj) and approximate 
-Vi « (Kl) ^ KlJA in and « {K^) « KlJ2 
in the elastic constants, where the average (•) was taken 

with respect to a circular Brillouin zone. The optimal 

Z^lfingth 03 is chosen such that Ed is minimal for a fixed 
amplitude ~ a corresponding to a typical defect elon- 
gation. * 

In the following, we treat first the more isotropic 
square vortex crystal YBCO (a = y^<j>o/B). From 
cqq and C44 for YBCO, the optimal length is given by 
03 = 4aAab/AcC(l — by/^. When comparmg the melt- 
— ipg criterium of the defect model in Eq. (|in|l with the 
^^indemann criterium obtained by equating the parame- 
ter Q to a fixed number, we get coincidgnce when tak- 
ing ii|djp account that the integrand in ^ and in lee of 
Eq. ^ receive their main contribution from the region 
~ Kbz/V2. We can put k^ fa in this region 
a^ceedfj? Ciii ~ 4/7r. Using this approximation 



SOrl esultingj 

in Eqs. and we can perform the integrals nunMg-- 
ically. Then we obtain from the melting condition llllu)) 
precisely ilae Lindemann criterium with the Lindcmann 
number fixed to be 
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Denoting by the spacing between the Cu02 double 
layers we obtain for the entropy jump per double layer 
and vortex 



A^z « kBT,^{dldT^){a,la^) ln[Zf| -\ (14) \ W2 
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Inserting ^ and (j^ we obtain 
ksTmas d 
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Finally, we make use of the Clausius-Clapeyron equation 
ASias/vQs — —{dHm/dT)AB/4TT, relating the jump of 
the entropy density ASia^/vQs to the jump of the mag- 
netic induction AB across the melting transition. Here 
Hm is the external magnetic field on the melting 1^ 
FromJjbe the Clausius-Clapeyron equation and Eqs. 
and nHl)|l ■ we obtain the following relations near Tc 
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FIG. 2: Entropy jump per layer per vortex AS; (first row) 
arid jump of magnetic induction field AB (second row) at the 



These results agree with the scaling results in Ref. I(TT>J7" 
where multiplicative prefactors are determined by fits to 
experimental curves (there is only a slight discrepancy 
because we use a different temperature dependence of 
the penetration depth). 

Next, we calculate the corresponding expressions in 
the case of the more layered crystal BSCCO (a ~ 
(2i/2/3i/'')v^<^o/S). First, we have to determine the dis- 
location length 03 in this case. For dislocation moves 
we have {u^y^^ ~ I/A'bz-* This mcgps that we can 
neglect the last two terms of C44 in nllz(l . coming from 
the self-energy of the vortex line, for the determination 
of the length of the dislocation. Remembering this we 
obtain by a similar procedure as for YBCO the dislo- 
cation length 03 fa 4:a\/2 Xab / Xcy/n . From this we find 
03065/0^044 <C 1, resulting in Iqq « 0. By taking into 
account that i?7r'^A^j/320o ^ 1 on the melting line we 
obtain that c, 
last term in 
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l^gi^ lting tr£|iijj^^ip|r^^ yhe experirtigf^tal values for YBCO are 
from Ref. (!l2|To^^^§jjf^d Ref. JiS for AB by squid experi- 
ments and Ref. py by torque measurCit^cm^^The experimi 
tal values for BSCCO are from Ref. Eu (squid) and Ref. 
(torque). 
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With the help of our melting criterion l|lin|l and 
the Clausius-Clapeyron equation (w(;_,use^diJ„i/dT « 
dB„JdT based on B,„(r) > HdiT) Iwfjf^e obtain for 
BSCCO: 
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By Hging these values we obtain by numerical integration 
of ||^ 
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Parameter ¥ahi^^ for optimal doped YBCO (BSCCO) 
are given by ff6f^afc(0) « 1186A (Aafc(O) « 2300A), 
Ca6(0) « 15A" (Ca6(0) « 30 A), the Cu02 double layer 
spacing = 12A (a^ = UA), = 92.7K (T^ = 90is:) 
13§r[ d the anisotropy parameter 7 = Xc/Xab « 5 (7 
200). We now calculate numerically the melting curves, 
the associated Lindemann parameter cl, the entropy and 
the magnetic induction jumps ASi and AB from the in- 
tersection crijterium ofthc full low- and high-temperature 
expressions (0' and without further approximations. 
The results are shown in Fjg. 1 anctPig. 2. Our ap- 
proximate analytic results ipiti|l and ipi9|l turn out to give 
practically the same curves. For comparisons, we show 

bQ|fciM^^wfegiJffi£Wital curves fw^^BC^^d^ 
14Rj:f. IXITliar^ and tor BSCCO of Ref. \m7Slt^^ — 



i^ab\ 



cj,laHog{l/cl) A, 



The first term comes from the integration region \k^\ < 
1/A|^j^, the second from the region 1/Xab < 1^31 < tt/os 
in 



temperature dependence of the Lindemann parameter for 
BSGCp seen in Fig.l comes mainly from the second term 
in IjllSII . which is the largest near T^,- Our curves are in 
reasonable agreement with the experimental data except 
at the end point of the melting line at low temperatures 
and near T « T^. That our vortex lattice model is not 
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a good approximation in these ranges is well known. At Near T k. T,. our model does not include the increase 
the low temperature side the discrepancy.^orncs^mainly of the entropy by thermal creation Qf,yortices beside the 
from Josephson decoupling of the layers iTlf most pro- magnetical ones forming the lattice l^f- For YB(^O ^glso. 
nounccd for the large anisotropic BSCCO supcrcondyc-^^^j^ order parameter fluctuations become important |21|r 
tor which results then also on large pinning effects \T2f . 
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We set up a melting model for vortex lattices in high-temperature superconductors based on the 
continuum elasticity theory. The model is Gaussian and includes defect fluctuations by means of 
a discrete-valued vortex gauge field. We derive the melting temperature of the lattice and predict 
the size of the Lindemann number. Our result agrees well with experiments for YBa2Cu307_5, 
and with modifications also for Bi2Sr2CaCu2 08 . We calculate the jumps in the entropy and the 
magnetic induction at the melting transition. 

PACS numbers: 74.25.Qt, 74.72.-h 



I. INTRODUCTION 

The magnetic flux lattices in high-temperature super- 
conductors can undergo a melting transition as was first 
suggested by Nelson in 1988 previously Brezin et al. 
[Si] had calculated a first-order liquid to solid phase tran- 
sition by renormalization group methods Q . Since then 
detailed properties of this transition have been studied 
in various theoretical and experimental papers. 

Most prominent are computer simulations methods 
based on the Langevin equation ,4] for the dynamics 
of the vortices or Monte Carlo simulations on the XY 
model coupled to an external magnetic field. Ana- 
lytic approaches are based mainly on the Ginzburg Lan- 
dau model ;G] or elastic models of the vortex lattice. 
The simplest estimates for the transition temperature 
in the vortex lattice came from an adaption of the fa- 
mous Lindemann criterion of three-dimensional ordinary 
crystals 0. Due to Houghton et al. and Brandt fsl, 
the criterium states that the vortex lattice undergoes a 
melting transition once the mean thermal displacement 
(m^)^/^ reaches a certain fraction of the lattice spacing 
a ~ ($0/^)^^^, where $0 is the flux quantum and B 
the magnetic induction. The ratio cl = {v?)^^"^ ja is the 
characteristic Lindemann number, which should be inde- 
pendent of B. Its value is not predicted by Lindemann's 
criterion. It must be extracted from experiments, and is 
usually found to lie in the range cl ~ 0.1 — 0.3. The 
most prominent examples of high-temperature super- 
conductors which exhibit vortex lattice melting are the 
anisotropic compound YBa2Cu307_5 (YBCO), and the 
strongly layered compound Bi2Sr2CaCu208 (BSCCO). 
Decoration experiments on BSCCO ^ show the forma- 
tion of a triangular vortex lattice, neutron scattering on 
YBCO of a tilted square lattice of vortices ^3 close to 
the melting region, the latter being favored by the d-wave 
symmetry of the order parameter and the anisotropy of 
the crystal luj . An explicit calculation of the Lindemann 
number cl = (u^)^/^/a for YBCO can be found in Ref. 
i and for BSCCO in Ref. [H 

In this paper, we present a theory which is capable of 
specifying the size of the Lindemann number cl. Ba- 
sis is a simple Gaussian model which takes into account 
both lattice elasticity and defect degrees of freedom in 



the simplest possible way. The relevance of defect fluc- 
tuations for the understanding of melting transitions is 
well-known. For ordinary crystals, this is textbook ma- 
terial In the context of vortex melting it was em- 
phasized in Ref. For ordinary crystals, the size of 
the Lindemann number has been calculated successfully 
by means of Gaussian lattice models with elastic and de- 
fect fluctuations, which clearly display flrst-order melting 
transitions in three dimensions and both first-order or se- 
quence of continuous transitions in two dimensions. 

An important quality of these models is that in the 
flrst-order case, where fluctuations are small, they lead 
to a simple universal melting formula determing the melt- 
ing point in terms of the elastic constants. The univer- 
sal result is found from a lowest-order approximation, in 
which one identifies the melting point with the intersec- 
tion of the high-temperature expansion of the free en- 
ergy density, dominated by defect fluctuations with the 
low-temperature expansion dominated by elastic fluctu- 
ations. The resulting universal formula for the melt- 
ing temperature determines also the size of the Linde- 
mann number. Recently, the results of Ref. for square 
crystals were successfully extended to face-centered and 
body-centered cubic lattices in three dimenisons ^] and 
also to two-dimensional triangular lattices A sim- 

ilar intersection criterium was also used for the melting 
point of vortex lattices in the Abrikosov approximation of 
the Ginzburg-Landau model useful for YBCO. Here 
we shall apply our defect model to calculate the melt- 
ing curves, the entropy jumps, and magnetic induction 
jumps of the vortex lattices in YBCO and BSCCO. We 
do not discuss in this work effects on the vortex lattice 
from other sources than defect fluctuations which can 
give rise to tricritical points and glass transitions T^,'!?!. 
Theoretical work on this subject can be found in Refs. 

There is in principle no problem of adding pin- 
ning in our formalism. For simplicities, we shall confine 
our discussion to the defect mechanism of melting. 

The melting criterium will be derived in Section II. 
The calculation of the melting temperatures, the entropy 
jumps, and the jumps of magnetic induction jumps for 
YBCO and BSCCO is carried out in Section III. 
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II. MELTING CRITERIUM 

Due to the large penetration depth Xat in the layers 
in comparison to a we have to take into account the full 
non-local elasticity constants when integrating over the 
Fourier space, as emphasized by Brandt in Ref. |^ For 
our Gaussian model, the partition function of the vortex 
lattice can be split into Z = ZgZfi where Zq is the parti- 
tion function of the rigid lattice and Zfi is thermally fluc- 
tuating part calculated via the elastic Hamiltonian plus 
defects. Due to the translational invariance of the vor- 
tex lattice in the direction of the vortices, which we shall 
take to be the 2;-axis, we may simply extend the models 
on square Q and triangular lattices by a third di- 
mension along the 0-axis, which we artificially discretize 
to have a lattice spacing 03, whose value will be fixed 
later. The elastic energy is 

Eel = ^ J2^y^U^){cll - 2C66)(V,U,) (1) 

X 

+ + VjmO C66 (ViUj + VjU.i) -I- {dzUi) C44 {dzUi). 

The subscripts i,j values 1,2 and l,m,n have values 
1, . . . , 3. The vectors Ui(x) are given by the transverse 
displacements of the line elements of the vortex lines with 
coordinate x. We have suppressed the spatial arguments 



of the elasticity parameters, which are really functional 
matrices Cy (x, x') = Cij(x — x'). Their precise forms 
were calculated by Brandt . The volume of the fun- 
damental cell V is equal to a^a^ (square) or 0^03-^/3/2 
(triangular). For a square lattice, the lattice derivate 
Vi in (HJ are given by Vi/(x) = [/(x + ae^) - /(x)]/a 
and V3/(x) = [/(x + 0363) - /(x)]/a3. For a triangu- 
lar lattice, the xy-pa,rt of the lattice has the link vec- 
tors ±ae(m) with 6(13) = (cos 27r/6, ± sin27r/6, 0) and 
6(2) = (~1: Oj 0). The lattice derivatives around a plaque- 
tte are defined by V(i)/(x) = [/(x -I- ae(i)) - /(x)] /a, 
V(2)/(x) - [/(x)-/(x-ae(2))]/a, V(3)/(x) = 
[/(x — ae(2)) — /(x -|- ae(i))] /a. From these we define 
discrete cartesian lattice derivatives used in the Hamil- 
tonian (|TJ V,/(x) = (2/3)e(,),V(,)/(x) and V3/(x) - 
[/(x -1-0363) — /(x)]/a3 transforming like the continuum 
derivative with respect to the symmetry group of the 
lattice 01 . Therefore, the Hamiltonian Q has the full 
symmetry of the triangular lattice and the correct con- 
tinuum elastic energy for zero lattice spacing. 

Within the elastic approximation the displacement 
fields are restricted to values within the fundamental cell. 
In order to contain also defect degrees of freedoms one 
has to insert into ^ integer valued defect gauge fields 
0. In the canonical stress representation, the partition 
function containing these fields becomes 



Zfi = det 



C66 



1/2 

det 



2tt(3 



4(cii - cge) 



5/2 



n 



n 



i<m 
Cii - 2C66 



dCTi, 



n E 

m nm(x) — — 00 , 





" du' 








exp 1 




J -00 ^ 





y- 



4(cii - cee) 



(E^+E'^'3^a.3 



(2) 



r 



The parameter (3 is given hy P — v CQQ/ksT {271)"^ 



give 

represent stress fields ^ . The matrix Dij (x) in Eq. (0) is 
a discrete- valued local defect matrix composed of integer- 
valued defect gauge fields m, 712,713 for square and 
triangular vortex lattices 0| as follows: 



D'7 



ni 71-3 
"3 "■2 



D 



A 



(3) 

-2^^l - "2 / 



The vortex gauge fields specify the VoUerra surfaces 
in units of the Burgers vectors. By summing over all 
51^1, 2, 3(x), the partition function Za includes all defect 
fluctuations, dislocations as well as disclinations. There 
is a constraint for a vortex lattice which does not exist 
for ordinary three-dimensional lattices. Dislocations in 
the vortex lattice can be both screw or edge type, but 



in either case the defect lines are confined in the plane 
spanned by their Burger's vector and the magnetic field 
13J, l2j|. The reason is that the flux lines in a vortex 
lattice cannot be broken. This results in the constraint 
£'11 = —D22 on the defect fields. 

We now calculate the low-temperature expansion of the 
partition function Zfi to lowest order, which includes only 
the Um = 0-term. By carrying out the integration over 
the displacement fields Mi(x) in (0) we obtain, as in 0, 
ITsl l. the leading term in the low-temperature expansion 
of the free energy 

1 

e 



7^- 



2N 



a ) det [(27r/3)c44/c66] 



hi 



where 



1 1 
2V^z 



- -— / d kdk-i log 



BZ 



C44 



.K*K, 



(5) 



(6) 
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Here Km is the eigenvalue of iVm- The momentum inte- 
grations in © run over the BriouUin zone of the vortex 
lattice whose volume is Vbz = (27r)'^/u, as indicated by 
the subscript BZ. 

Next we calculate the high-temperature expansion 
to lowest order. By carrying out the integration 
over the displacement fields Ui(x) in Q and further by 
summing over the defect fields rtm under the condition 
D\\ — — 1?22 mentioned above, it turns out that the stress 
fields a\2 and cr„ = a\\ — 022 can have only integer num- 
bers. Doing the integrals over the stress fields (7^3 and 
Ug = ail -\- C22 we obain, in the lowest order high tem- 
perature limit, corresponding to (T12 = and (t„ = 0: 



a 



N 



-Nh 



2^ det [(2^/3)2c44/C66] 



with 



1 1 
2Vb^ 



d kdkz log 



BZ 



Cll - C66 



C44 KXK^ 



(7) 



(8) 



The constant C has the values = 1 for the square 
vortex lattice and Ca = \/3 for the triangular one. 

In the low-tempcraturc expansion representing the 
solid phase, defect field configurations n„i ^ correspond 
to dislocations and disclinations giving finite temperature 
corrections to the free energy —\og{Z)/kBT. These cor- 
rections are exponentially small with an exponent pro- 
portinal to —f3 0,01 • In contrast to this corrections to 
the high temperature expansion in the fiuid phase cor- 
responding to stress configurations CT12 ^ and (Tu ^ 
of integer values result in temperature corrections to the 
free energy which are also exponentially small with an 
exponent proportional to — 1//3. The structure of the 
high- and low-temperature corrections to the partition 
function is extensively discussed in Refs. 0, E for ordi- 
nary crystals, and can be easily transferred to our case 
of vortex lattices. It was shown in Refs. 0,^^for the two 
dimensional square and triangular as well as the three di- 
mensional square crystal that the exponentially vanishing 
higher order corrections to the low- and high-temperature 
expansion of the free energy are negligible in the determi- 
nation of the transition temperature. This is particularly 
true for the three dimensional crystal (see p. 1082 in 0) 
which we take as a justification to restrict our calculation 
to lowest order in this paper. 

From the partition function (0) with no defects {n„i — 
0) we obtain for the Lindemann number = (?i^)^/^/a 
the momentum integral 



knT 



a^v Vbz Jbz 



d-kdk-i 



^^^,t^K*K,+alKIK, 
(9) 

This can be simplified by taking into account that cii is 
much larger than cge, C44 in the relevant regime |2ll. l23j. 
As announced, we find the melting temperature from the 



intersection of low-and high-temperature expansions, ob- 
tained by equating Z'^^'^ = By taking into ac- 
count det[a|V3V3] = 1 we obtain h,lii ^ Iqq, and fur- 
ther that the i = 1-term in © is much smaller than the 
i — 6-term. In the following analytic discussion (but not 
in the numerical plots) we neglect these small contribu- 
tions. The temperature of melting is then given by the 
simple formula 



knT 



1 



V det 



l/N 



■C 



(10) 



where det[c66] is the determinant of the N x N functional 
matrix cgg . The elastic moduli C44 and cge at low reduced 
magnetic fields b = B/Hc2 < 0.25 can be taken from 
Brandt's paper [23| 



C66 



(11) 



C44 = 



ab 3 



47r(l + A2P + A2,fc32) 32^2^2 i + A^i^BZ + A^^fcl 

^32n^Xt,ki l + Xl,kl/{l + 2Xljm ^ > 

where Ac is the penetration depth in the xy-plane, C = 1, 
and Kbz is the boundary of the circular Brillouin zone 
^BZ — 47ri?/0o- At high fields (b > 0.5), cee is altered 
by a factor ( « 0.71(1 — b), and the penetration depths 
in C66, C44 are replaced by A^ = A-^/(l — 6), where A de- 
notes either Xab or Ac. In addition, the last two terms in 
C44 are replaced by B^o/IStt^A^. For YBCO we have 111 

A(r) = A(o)[i-(r/rc)]-V3, ^(t) = c(o)[i-(T/rc)]-i/2 

for BSCCO A(T) = A(0)[1 - {T/T^)'^]-^^^, ^{T) = 
C(0)[1 - (T/rc)4]-i/2 . When calculating C44 in ^ we 
have used a momentum cutoff in the two-vortex inter- 
action potential k < 2/(m^)^/^, and not the inverse of 
the correlation length 1/^ as in Ref. |^ The cutoff is 
due to thermal softening and becomes relevant for 
> 1, or equivalently for clV2^^/Hc2{.T)/B > 
1, which is fulfilled in the melting regime of BSCCO, but 
not for YBCO. 

It remains to determine the effective lattice spacing 
03 along the vortex lines. An elementary defect in 
the vortex lattice (arising for example from a crossing 
of two vortex strings) takes place over a typical length 
scale in the z-direction determined from the condition 
that the sum of elastic displacement energy and the en- 
ergy required to stretch the line against the line ten- 
sion is minimal. It is the elastic energy of this small- 
est defect that has to be taken into account for in our 
model. The energy of an ensemble of dislocations is de- 
termined by the interplay of elastic energy of small dis- 
placements and integer-valued defect fields. The rele- 
vant part of the free energy is given by the discretized 
free energy — log(Zo-^fl) fcs?^ (El) in which 03 is equal 
to the above length scale in the z-direction. To deter- 
mine this, we insert the variational ansatz for the trans- 
verse displacement field Ui = J^^iAq exp[— 2|z|/a3] into 
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the continuum version (in z direction) of the elastic en- 
ergy £3) and approximate — ~ ~ ^iz/^ i^i -^oi 
and (-^^) ~ -^bz/^ in the elastic constants, where 
the average (...) was taken with respect to a circular 
Brillouin zone. The optimal length scale 03 is chosen 
such that £^01 is minimal for a fixed amplitude Aq ~ a 
corresponding to a typical defect elongation. 



III. APPLICATION TO YBCO AND BSCCO 

In the following, we treat first the more isotropic 
square vortex crystal YBCO (a — ^Jcf)Q/B). From cee 
and C44 for YBCO, the optimal length scale is given by 
as = 4aAafc/AcC(l — 6)^^^. When comparing the melting 
criterium of the defect model in Eq. H1U|) with the Linde- 
mann criterium obtained by equating the parameter ^ 
to a universal number, we obtain identical results when 
taking into account that the integrand in 1^1 and in /ge 
of Eq. 10 receive their main contribution from the region 
k « \J (fc^) « Kbz/V^- We can approximate ^3 « in 
this region, resulting in a\cQ(, / o?' ~ A:/ it. With the 
same approximation in Eqs. lO and lO, we can perform 
the integrals numerically. Then we obtain from the melt- 
ing condition precisely the Lindemann criterium in 
which the Lindemann number (|5J| is predicted to be 



BJ-m 



C44(^,0)C66(^,0) 



1/2 



(0.18)2. (13) 



Denoting the spacing between the Cu02 double layers by 
as we obtain for the entropy jump per double layer and 
vortex 



« kBTra{d/dT„,){as/a^) ln[Z^'^°°/Zj^°]. (14) 



Inserting lO and Q, this becomes 



ksTrnas d 



03 



dT„ 



In- 



.^44 1 ^2 



■,0)C66('^' 



72 



,0) 



(15) 

Finally, we make use of the Clausius-Clapeyron equation 
which relates the jump of the entropy density across the 
melting transition to the jump of the magnetic induc- 
tion by ASia^/vas = -{dHm/dT)AB /An. Here iJ^ is 
the external magnetic field on the melting line. Combin- 
ing the Clausius-Clapeyron equation with Eqs. H12I) and 
(|15|l we obtain, with the abbreviation Tm = Tm/Tc, the 
following relations near Tci 



AB 



12 c (i-r/rj4/3e4 05 

le^TT^ {ksTrxi^mm 

VCas Ac fc_B _ 2/7 _ 
6a Aq6(1 

VCtT Ac 



a. c 



2 ja 



s ^Ly'o 



T„.) 103 Te(l - T„ 



2.5 (1 



\2/3 2 



^/^A2,(G)' 
(16) 



These results agree with the general scaling results in 
Ref. [23| , with the advantage that here the prefactors are 
predicted whereas those in had to be determined by 
fits to experimental curves (there is only a slight discrep- 
ancy because we use a different temperature dependence 
of the penetration depth). 

Next, we calculate the corresponding expressions in 
the case of the more layered crystal BSCCO (a = 

(21/2/3^/"*) v/0o7s)- First, we have to determine the 
dislocation length scale in this case. For dislocation 
moves wc have (u^)^/^ ~ 1/Kbz- This means that we can 
neglect the last two terms of C44 in H12|l , coming from the 
self-energy of the vortex line, when determining as. Re- 
membering this we obtain by a similar procedure as for 
YBCO the dislocation length scale « 4a-\/2 Aab/AcV^- 
From this we find a|c66/a2c44 ^ 1, resulting in Iqq « 0. 
By taking into account that 371^X^1^/32(1)0 < 1 on the 
melting line we obtain that C44(A;,fc3) for jfesj < tt/as is 
dominated by the last term in 112|l . Then we obtain 



C44(A;, fcs) 



32n^Xl,{l + Xl,K, 



for fcs < , (17) 



a . ^ ^ i?0oln(l + 2i?Ay0oci) , , 1 

" — 32^^5Api — » a:;- 



By using these values we obtain by numerical integration 
of® 



ksTm ■ 0.36 



^ ksTmas ■ 1.60 

(^,0)C44(^,0) a^C44(^,l/a3) 



V2 



(18) 

^ ksTrnXl^ ■ 138 /-I , \2 1^2 I ksTjnXl^Xl ■ 137 Xab 

<^la y^ + ^at^BZ+ ^2aHogil/cl) A.- 

The first term comes from the integration region jfcsj < 
l/Xab, the second from the region l/Xab < l^sl < t^/o-3 
in ®. 

From our melting criterion (|10|l and the Clausius- 
Clapeyron equation (where dHm/dT w dBm/dT due to 
Bm{T) > Hci{T) 13), we obtain for BSCCO 



ASi 



AB 



1 1(1- (T/Tc)4) 



192^3^7 AL(0)A2(0) {kBTf 
y/Trasks Ac 1 + Sr^ 2.9 0^(1 -t- 3r, 



4 \a2 



4y2a Xab 1-t4 104 T,„A2,(0) ' 

4 



7^'/' ^c. rj. ^ 1-8 (1 - r„^m 



2V2a Xab 



10^ A2,(0) 



(19) 



Parameter values for optimal doped YBCO (BSCCO) 
are given by ^ Aab(O) « 1186A (Aa6(0) w 2300A), 
£,ab{0) ~ 15 A {(,ab{0) « 30A), the CuOz double layer 
spacing = 12A (a^ = 14A), Tc = 92.7 K (Tc = 90K) 
and the anisotropy parameter 7 = Xc/Xab ~ 5 (7 « 200). 

We now calculate numerically the melting curves, the 
associated Lindemann parameter cl, the entropy and the 
magnetic induction jumps ASi and AB from the inter- 
section criterium of the full low- and high-temperature 
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FIG. 1: Melting curve B = Bm{T) for YBCO and BSCCO. 
The exp erimental values are from Ref. for YBCO and 
Ref. [23 for BSCCO. The numbers on the theoretical melting 
curves are the Lindemann numbers cl calculated from J^J. 
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FIG. 2: Entropy jump per double layer per vortex ASi (first 
row) and jump of magnetic induction field AB (second row) 
at the melting transition. The experimental values for YBCO 
are from Ref. HI for AS';, Ref.lp for AB by SQUID ex- 
periments (SQUID), and Ref. |33 by torque measurements 
(Torque) . The experimental values for BSCCO are from Ref. 
I31] by SQUID measurements (SQUID) of the magnetic field 
and Ref. |23by microscopic Hall sensors (Hall). 

expressions (O and (|7|l without further approximation. 
To accomplish this, we use the elastic constant cn given 
by Brandt in Ref. The intersection criterium 

of the low and high temperature expansion of the 
partition function is then at least in the case of BSCCO 
a complicated integral equation via the dependence of 



C44 on the Lindemann parameter cl- One can solve this 
integral equation by numerical methods. The results are 
shown in Fig. 1 and Fig. 2. 

IV. DISCUSSION 



Our approximate analytic results Hlt)|) and H19|) for 
YBCO and BSCCO turn out to give practically the same 
curves. For comparisons, we show in bot h fig ures the ex- 
perimental curves for YBCO of Ref. El El 113 and for 
BSCCO of Ref. El EH We see in Fig. 1 and Eq. ^ 
that the Lindemann number is independent of the mag- 
netic field for YBCO for the entire temperature range. 
For BSCCO we obtain a magnetic field dependence of 
the Lindemann parameter near Tc- This comes mainly 
from the second term in H18(l . having its largest variation 
near Tc- The theoretical melting curves in Fig. 1 and the 
entropy and the magnetic induction jumps in Fig. 2 are 
in reasonable agreement with the experimental curves, 
except at the end point of the melting line at low tem- 
peratures and near T k T^. It could have been antici- 
pated that our vortex lattice model is not a good approx- 
imation in these regimes. At low temperature, the dis- 
crepancy comes mainly from Josephson decoupling of the 
layers |32| |. most pronounced for the strongly anisotropic 
BSCCO superconductor, which leads also to large pin- 
ning effects [s^. We think that this is also the reason 
for the difference in the curves of Kadowki et al. in Ref. 
I31I and of Zeldov et al. in Ref. E3 shown in Fig. 2. Pin- 
ning has the largest influence on the form of the melting 
curve at low temperatures |34l | resulting in a decrease of 
AS"; and AJ5 j23| in the limit of low temperatures shown 
by the curves of Zeldov et al.. Near T k, Tc, our model 
does not include the increase of the entropy by thermal 
creation of vortices in addition to the ones caused by the 
external magnetic field which form the lattice [33]. For 
YBCO, also order parameter fluctuations become impor- 
tant |36| which are ignored here. 

Summarizing, we obtain within our formalism the 
melting curve, the entropy and the magnetic jump within 
a single theory. This is not possible by the widely used 
Lindemann rule. Our curves agree well with the experi- 
mental curves for YBCO and BSCCO except at the end- 
points. We gave analytical approximations for the var- 
ious quantities. Our theory is a direct generalization of 
the elastic theory of vortex displacements by taking into 
account the defect degrees of freedom of the lattice which 
makes it possible to obtain the physics of the liquid phase 
of the vortex lattice beside the sohd phase. 
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and with modifications also for Bi2Sr2CaCu208 . We calculate the jumps in the entropy and the 
magnetic induction at the melting transition. 
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I. INTRODUCTION 

The magnetic flux lattices in high-temperature super- 
conductors can undergo a melting transition as was first 
suggested by Nelson in 1988 Previously, Brezin et 

al. Q had calculated a first-order liquid to solid phase 
transition by renormalization group methods Since 
then detailed properties of this transition have emerged 
from various theoretical and experimental papers. 

Most prominent are computer simulations of the 
Langevin equation Q for the dynamics of the vortices, 
or Monte Carlo simulations of XY-type model Q cou- 
pled to an external magnetic field. Analytic approaches 
are based mainly on the Ginzburg-Landau model |6|, or 
on elastic models of the vortex lattice. The simplest esti- 
mates for the transition temperature in the vortex lattice 
came from an adaption of the famous Lindemann crite- 
rion of three-dimensional ordinary crystals ■ In the for- 
mulation of Houghton et al. and Brandt , the criterium 
states that the vortex lattice undergoes a melting transi- 
tion once the mean thermal displacement (m^)^/^ reaches 
a certain fraction of the lattice spacing a « (<i>o/i3)^/^, 
where $0 is the flux quantum and B the magnetic induc- 



tion. The ratio cl 



,2\l/2 



/a is the characteristic Lin- 



demann number, which should be independent of B. Its 
value is not predicted by Lindemann's criterion. It must 
be extracted from experiments, and is usually found to 
lie in the range cl ~ 0.1 — 0.3. 

The most prominent examples of high-temperature 
superconductors which exhibit vortex lattice melting 
are the anisotropic compound YBa2Cu307_5 (YBCO), 
and the strongly layered compound Bi2Sr2CaCu2 08 
(BSCCO). Decoration experiments on BSCCO % show 
the formation of a triangular vortex lattice, neutron scat- 
tering on YBCO of a tilted square lattice of vortices 
close to the melting region, the latter being favored by 
the d-wave symmetry of the order parameter and the 
anisotropy of the crystal explicit calculation of 



the Lindemann number cl 



,2x1/2 



/a for YBCO can 



be found in Ref. |^ and for BSCCO in Ref. [H 

In this paper, we present a theory which is capable of 
specifying the size of the Lindemann number cl, and pre- 



dicting corrections to the criterium. Our theroy is based 
on a simple Gaussian model which takes into account 
both lattice elasticity and defect degrees of freedom in 
the simplest possible way. The relevance of defect fluc- 
tuations for the understanding of melting transitions is 
well-known. For ordinary crystals, this is textbook ma- 
terial In the context of vortex melting it was em- 
phasized in Ref. Foi' ordinary crystals, the size of 
the Lindemann number has been calculated successfully 
by means of Gaussian lattice models with elastic and de- 
fect fluctuations, which clearly display flrst-order melting 
transitions in three dimensions and both first-order or se- 
quence of continuous transitions in two dimensions. 

An important quality of these models is that in the 
first-order case, where fluctuations are small, they lead 
to a simple universal melting formula determing the melt- 
ing point in terms of the elastic constants. The univer- 
sal result is found from a lowest-order approximation, in 
which one identifies the melting point with the intersec- 
tion of the high-temperature expansion of the free en- 
ergy density, dominated by defect fluctuations with the 
low-temperature expansion dominated by elastic fluctu- 
ations. The resulting universal formula for the melt- 
ing temperature determines also the size of the Linde- 
mann number. Recently, the results of Ref. for square 
crystals were successfully extended to face-centered and 
body-centered cubic lattices in three dimenisons and 
also to two-dimensional triangular lattices A sim- 

ilar intersection criterium was also used for the melting 
point of vortex lattices in the Abrikosov approximation of 
the Ginzburg-Landau model 0] useful for YBCO. Here 
we shall apply our defect model to calculate the melt- 
ing curves, the entropy jumps, and magnetic induction 
jumps of the vortex lattices in YBCO and BSCCO. We 
do not discuss in this work effects on the vortex lattice 
from other sources than defect fluctuations which can 
give rise to tricritical points and glass transitions [l^llTj . 
Theoretical work on this subject can be found in Refs. 
0,0,123 There is in principle no problem of adding pin- 
ning in our formalism. For simplicities, we shall confine 
our discussion to the defect mechanism of melting. 

The melting criterium will be derived in Section II. 
The calculation of the melting temperatures, the entropy 



2 



jumps, and the jumps of magnetic induction for YBCO 
and BSCCO is carried out in Section III. 



II. MELTING CRITERIUM 

Due to the large penetration depth Xat in the layers 
in comparison to a we have to take into account the full 
non-local elasticity constants when integrating over the 
Fourier space, as emphasized by Brandt in Ref. For 
our Gaussian model, the partition function of the vortex 
lattice can be split into Z ^ ZaZa where Zq is the parti- 
tion function of the rigid lattice and Zfi is thermally fluc- 
tuating part calculated via the elastic Hamiltonian plus 
defects. Due to the translational invariance of the vor- 
tex lattice in the direction of the vortices, which we shall 
take to be the 2;-axis, we may simply extend the models 
on square [rj and triangular lattices [Tsj by a third di- 
mension along the z-axis, which we artificially discretize 
to have a lattice spacing 03, whose value will be fixed 
later. The elastic energy is 

^(ViUi)(cii - 2c66)(ViUi) (1) 

X 

+ ^(ViMj + WjUi) C66 (ViUj + VjUi) -f (VaUi) C44 (VaU;). 

The subscripts i,j have values 1, 2 and l,m,n have values 
1, . . . ,3. The vectors Ui(x) are given by the transverse 
displacements of the line elements of the vortex lines with 
coordinate x. We have suppressed the spatial arguments 
of the elasticity parameters, which are really functional 
matrices Cy (x, x') = Cij(x — x'). Their precise forms 
were calculated by Brandt . The volume of the fun- 
damental cell V is equal to 0^03 (square) or a^a3-\/3/2 
(triangular). For a square lattice, the lattice derivate 
Vi in (UJ are given by Vi/(x) = [/(x -I- ae^) - /(x)]/a 
and V3/(x) = [/(x -|- 0363) - /(x)]/a3. For a triangu- 
lar lattice, the xy-part of the lattice has the link vec- 
tors ±ae(m) with e(i 3) = (cos 27r/6, ± sin27r/6, 0) and 
G(2) = ("Ij 0, 0). The lattice derivatives around a plaque- 
tte are defined by V(i)/(x) = [/(x + ae(i)) — /(x)] /a, 
V(2)/(x) = [/(x)-/(x-ae(2))]/a, V(3)/(x) = 
[/(x — ae(2)) — /(x -I- aej-i))] /a. From these we define 
discrete cartesian lattice derivatives used in the Hamil- 



tonian O V,/(x) = (2/3)e(,),V(,)/(x) and V3/(x) = 
[/(x -I- 0363) — /(x)]/a3 transforming like the continuum 
derivative with respect to the symmetry group of the 
lattice (T^. Therefore, the Hamiltonian ^ has the full 
symmetry of the triangular lattice and the correct con- 
tinuum elastic energy for zero lattice spacing. 

The quadratic approximation to the energy ^ is so 
far only appropriate for the the low-temperature clas- 
sical thermodynamic behaviour. It is possible to ex- 
tend the Hamiltonian at the quadratic level in such a 
way that the range of applicability extends beyond the 
melting transition. This is possible by the introduction 
of integer- valued defect gauge fields. We observe that 
the displacement fields in are restricted to values 
within the fundamental cell. The defect gauge fields en- 
ter to characterize the jumps of the displacement fields 
across the Volterra surface 0, ^3 ■ As usual for gauge 
fields we choose a minimal coupling to the lattice dis- 
placements. On account of the three lattice derivates 
per fundamental cell and further two dimensions for the 
displacements there are six independent integer-valued 
gauge fields per fundamental cell corresponding to the 
various defect configurations. One can eliminate two of 
them (we choose the defect fields corresponding to jumps 
in the z-direction) by relaxing the restriction of the dis- 
placement fields to the fundamental cell (elimination of 
gauge freedom) . See the discussion in Ref. E for square 
lattices. A similar consideration was also carried out in 
[T5I I for the two-dimensional triangular lattice where the 
elimination of the gauge degrees of freedom is more com- 
plicated due to the absence of the ^-direction. Finally, we 
can eliminate one more integer-valued defect field since 
the elastic energy in depends only on the displace- 
ments Ui via the strain tensor ViUj -I- VjUi jTj]. In sum- 
mary, only three independent integer-valued fields have 
to be included in Q to obtain the elastic energy of the 
vortex lattice including defects. By taking into account 
the above considerations one can then easily determine 
the partition function including defects for the square 
vortex lattice by using the considerations in Ref. and 
for the triangular ones by using Ref. Il5i 

By using a Hubbard-Stratonovich decoupling of the 
quadratic displacement terms in (Q, the stress represen- 
tation of the partition function becomes 



Zfi = det 



C66 



4(cii - cee) 



1/2 
det 



27r/3 



5/2 



n 



n/ da. 



n E 

m iini(x) — — 00 



E 4 + ^ E 4 - (E 4(cn -'cI) (E + E ^-^3 < 

i<j i i ^ ^ i i } ) 



The parameter (3 is given by /3 = v ceQ/kBT{27r)'^ . aij represent stress fields W\. The matrix (x) in Eq. Q is 



3 



a discrete- valued local defect matrix composed of integer- 
valued defect gauge fields ni,n2,n3 for square and 
triangular vortex lattices as follows: 



"3 n2 



2^1 — "2 



(3) 



(4) 



The vortex gauge fields specify the Volterra surfaces 
in units of the Burgers vectors. By summing over all 
?^i,2,3(x), the partition function includes all defect 
fluctuations, dislocations as well as disclinations. There 
is a constraint for a vortex lattice which does not exist 
for ordinary three-dimensional lattices. Dislocations in 
the vortex lattice can be both screw or edge type, but 
in either case the defect lines are confined in the plane 
spanned by their Burger's vector and the magnetic field 
[l^ |2^ . The reason is that the flux lines in a vortex 
lattice cannot be broken. This results in the constraint 
Dii — —D22 on the defect fields. 

We now calculate the low-temperature expansion of the 
partition function Za to lowest order, which includes only 
the n„i = 0-term. By carrying out the integration over 
the displacement fields iii(x) in Q we obtain, as in 0, 
[l5l |. the leading term in the low-temperature expansion 
of the free energy 



0-3 



2N 



1 



a J det [(27r/3)c44/c66] 



£{1,6} ' 



(5) 



where 



11/"? 
= n T7~ / d kdkz log 
^ VBZ Jbz 



C44 



-K*K,+alKlK^ 



(6) 

Here Km is the eigenvalue of iWm-N is the number of 
vertices in the lattice. The momentum integrations in 
run over the BriouUin zone of the vortex lattice whose 
volume is Vbz — (27r)'^/w, as indicated by the subscript 
BZ. 

Next we calculate the high-temperature expansion 
to lowest order. By carrying out the integration 
over the displacement fields Wi(x) in (0) and further by 
summing over the defect fields rim under the condition 
Dii — —D22 mentioned above, it turns out that the stress 
fields (T12 and (t„ = an — G22 can have only integer num- 
bers. Doing the integrals over the stress fields aa and 
tjg = (Til + 1722 we obain, in the lowest order high tem- 
perature limit, corresponding to o\2 = and tT„ = 0: 



1 



with 



\a) 2N det [(2^/3)2c44/c66] 

cii - C66 K*Kj 



h^]--^ I d^kdks log 
2 Vbz Jbz 



1 + 



C44 k;k3 



(7) 



(8) 



The constant C has the values Cn = 1 for the square 
vortex lattice and Ca = a/S for the triangular one. 

In the low-temperature expansion representing the 
solid phase, defect field configurations rim 7^ correspond 
to dislocations and disclinations giving finite temperature 
corrections to the free energy —log{Z)/kBT. These cor- 
rections are exponentially small with an exponent pro- 
portinal to — /3 0, In contrast to this corrections to 
the high temperature expansion in the fluid phase cor- 
responding to stress configurations tTi2 9^ and CTu 7^ 
of integer values result in temperature corrections to the 
free energy which are also exponentially small with an 
exponent proportional to — The structure of the 
high- and low-temperature corrections to the partition 
function is extensively discussed in Refs. for ordi- 

nary crystals, and can be easily transferred to our case 
of vortex lattices. It was shown in Refs. 0,^1 for the two 
dimensional square and triangular as well as the three di- 
mensional square crystal that the exponentially vanishing 
higher order corrections to the low- and high-temperature 
expansion of the free energy are negligible in the determi- 
nation of the transition temperature. This is particularly 
true for the three dimensional crystal (see p. 1082 in 7]) 
which we take as a justification to restrict our calculation 
to lowest order in this paper. 

From the partition function ^ with no defects = 
0) we obtain for the Lindemann number cl = (u^)^^^ /a 
the momentum integral 



kuT 



a^v Vbz Jbz 



d kdk-i 



C44 



K*K. 



alK*K^ 



(9) 

This can be simplified by taking into account that Cn is 
much larger than cqq , C44 in the relevant regime |^ . 
As announced, we find the melting temperature from the 
intersection of low-and high-temperature expansions, ob- 
tained by equating Z^^^ — Z'^^°°. By taking into ac- 
count det[a|V3V3] = 1 we obtain h^lu <C ^66, and fur- 
ther that the i — 1-term in © is much smaller than the 
i = 6-term. In the following analytic discussion (but not 
in the numerical plots) we neglect these small contribu- 
tions. The temperature of melting is then given by the 
simple formula 



knT 



1 



V deti/^[c66: 



(10) 



where det[c66] is the determinant of the N x N functional 
matrix cge . The elastic moduli C44 and cge at low reduced 
magnetic fields b = B/Hc2 < 0.25 can be taken from 
Brandt's paper i21| 



C66 



C44 



(11) 



rln- 



,6'^3 



47r(l + A2p-HA2^fc2) 32^2^2-1 



\lKBZ + \l,kl 



+ 



i2n^K,kl 



In 



l + A^,fci/(H-2Ay(u)2) 



(12) 
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where Ac is the penetration depth in the xy-plane, C = li 
and K^z is the boundary of the circular Brillouin zone 
-^BZ = 4:TrB/(j)o. At high fields (6 > 0.5), cge is altered 
by a factor ^ « 0.71(1 — b), and the penetration depths 
in C66, C44 are replaced by = A^/(l — 6), where A de- 
notes either Xab or Ac. In addition, the last two terms in 
C44 are replaced by B(j)o/16'K'^Xl. For YBCO we have Q 
A(T) = A(0)[l-(T/rcri/3, ^(T) ^ ^(0)[l-(T/rc)]-i/2 

for BSCCO lH A(T) = A(0)[1 - (r/Tc)*]-!/^^ ^(T) = 
e(0)[l - (r/Tc)4]i/2/[i _ (T/Tc)2] . When calculat- 
ing C44 in H12|l we have used a momentum cutoff in 
the two- vortex interaction potential k < 2/(u^)^/^, and 
not the inverse of the correlation length 1/f as in Ref. 



|2U The cutoff is due to thermal softening |7|, and be- 
comes relevant for (u^)^/^/^ ^ 1, or equivalently for 
CLV2ny/Hc2{T)/B > 1, which is fulfilled in the melt- 
ing regime of BSCCO, but not for YBCO. 

It remains to determine the effective lattice spacing 
03 along the vortex lines. An elementary defect in 
the vortex lattice (arising for example from a crossing 
of two vortex strings) takes place over a typical length 
scale in the z-direction determined from the condition 
that the sum of elastic displacement energy and the en- 
ergy required to stretch the line against the line ten- 
sion is minimal. It is the elastic energy of this small- 
est defect that has to be taken into account for in our 
model. The energy of an ensemble of dislocations is de- 
termined by the interplay of elastic energy of small dis- 
placements and integer-valued defect fields. The rele- 
vant part of the free energy is given by the discretized 
free energy —log{ZQZfi) ksT ((SJ in which 03 is equal 
to the above length scale in the ^-direction. To deter- 
mine this, we insert the variational ansatz for the trans- 
verse displacement field Ui — Ji^iAg exp[— 2|z|/a3] into 
the continuum version (in z direction) of the elastic en- 
ergy Q and approximate — V2 ~ {K2) ~ K^^/i in E^i 
and w (-^^) ~ -^bz/^ the elastic constants, where 
the average (...) was taken with respect to a circular 
Brillouin zone. The optimal length scale 03 is chosen 
such that £^01 is minimal for a fixed amplitude Ao ~ a 
corresponding to a typical defect elongation. 



III. APPLICATION TO YBCO AND BSCCO 

In the following, we treat first the more isotropic 
square vortex crystal YBCO (a — y^^oAB)- From cqq 
and C44 for YBCO, the optimal length scale is given by 
03 = 4aAafc/AcC(l — by^^. When comparing the melting 
criterium of the defect model in Eq. H1U|) with the Linde- 
mann criterium obtained by equating the parameter 
to a universal number, we obtain identical results when 
taking into account that the integrand in ^ and in Zge 
of Eq. (O receive their main contribution from the region 
k « (fc2) « Kbz/V^- We can approximate A:3 « in 
this region, resulting in a^cee/ 0^044 « A/tt. With the 
same approximation in Eqs. and we can perform 
the integrals numerically. Then we obtain from the melt- 



ing condition (|10|l precisely the Lindemann criterium in 
which the Lindemann number ||^ is predicted to be 



kuT, 



BJ-m 



C44(^,0)C66(^,0) 



1/2 



(0.18)2. (13) 



Denoting the spacing between the Cu02 double layers by 
as we obtain for the entropy jump per double layer and 
vortex 



^Si « kBT^{d/dT^){as/a^)HZl-^/Zl 



(14) 



Inserting ||SJ) and {Tj), this becomes 



knT, 



as 



d_ 



In- 



ksTni/a? 



C44(^,0)C66(^,0) 



ATbz 



(15) 

Finally, we make use of the Clausius-Clapeyron equation 
which relates the jump of the entropy density across the 
melting transition to the jump of the magnetic induc- 
tion by ASia^/vag = -{dHm/dT)AB /in. Here Hm is 
the external magnetic field on the melting line. Combin- 
ing the Clausius-Clapeyron equation with Eqs. (|12() and 
(I15|l we obtain, with the abbreviation r,„ = Tm/Tc, the 
following relations near Tc: 



Bm{T) 
ASi 
AB 



12 c (i-r/rc)4/3ci03 

le^TT^ (fcBT)2A2,(0)A2(0) ' 
VCtts Ac kB 2.7 



a., c 



2 a2 



s ^lVQ 



6a A,b(l-r„0 103 rc(l - r„0i/3A2,(0) ' 

V^TT Ac , ^ 2.5 (l-T„)2/3c2 0o 



2a4>o Xab 



-knTrr 



102 



KM 



(16) 



These results agree with the general scaling results in 
Ref. [2^ , with the advantage that here the prefactors are 
predicted whereas those in |2^ had to be determined by 
fits to experimental curves (there is only a slight discrep- 
ancy because we use a different temperature dependence 
of the penetration depth). 

Next, we calculate the corresponding expressions in 
the case of the more layered crystal BSCCO (a = 
(2i/2/3i/'')v^(/)o/S). First, we have to determine the 
dislocation length scale 03 in this case. For dislocation 
moves we have (u^)^/^ ~ 1/ K^i- This means that we can 
neglect the last two terms of C44 in H12|l , coming from the 
self-energy of the vortex line, when determining 03. Re- 
membering this we obtain by a similar procedure as for 
YBCO the dislocation length scale 03 w 4a'\/2 Aab/Ac^/TT- 
From this we find a^CQQ/o?'C44 <C 1, resulting in ^ge ~ 0. 
By taking into account that BTi^X^^^/'i2(t)Q < 1 on the 
melting line we obtain that C44(A:,fc3) for \kj,\ < nja^i is 
dominated by the last term in H12|l . Then we obtain 

32TT-^Xl^k^ Xab 



5 



T[K] (YBCO) 




m 0.080.07 90 



T[K] (BSCCO) 



FIG. 1: Melting curve B = Bm{T) for YBCO and BSCCO. 
The exp erimental values are from Ref. for YBCO and 
Ref. [23 for BSCCO. The numbers on the theoretical melting 
curves are the Lindemann numbers cl calculated from J^J. 



By using these values we obtain by numerical integration 
of ® 



kRTm ■ 0.36 



kBTma-i ■ 1-60 



a^/C66(^,0)c44(^,0) a*C44(^,l/a3) 



kBTmKb ■ 138 

00 " 



1 + \luK- 



(18) 

ksTrnXl^l^X^ ■ 137 Xab 



The first term comes from the integration region \k^\ < 
1/Xab, the second from the region l/Xat < l^sl < Tr/aa 
in 10. 

From our melting criterion (|10|l and the Clausius- 
Clapeyron equation (where dHm/dT « dBm/dT due to 
B„i{T) > Hci{T) 13), we obtain for BSCCO 



AB 



1 



1 (1 - (T/T,)4)2 



192^3^7 A2^(0)A2(0) (fcBT)2' 
V^askB Xc 1 + 3r4 _ 2.9 0,(1 + 3x^)02 
4V2a Aafa 1 - t4 ^ 104 T™A2,(0) ' 



2V2a Xai 
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(19) 



Parameter values for optimal doped YBCO (BSCCO) 
are given by |2| Aab(O) « 1186A (Aa6(0) « 2300A), 
^a6(0) « I5A (Ca6(0) « 30A), the CuOa double layer 
spacing a, = 12A (a, = UA), = 92. (T^ = 90^) 
and the anisotropy parameter 7 = Xc/Xab ~ 5 (7 ~ 200). 

We now calculate numerically the melting curves, the 
associated Lindemann parameter cl , the entropy and the 
magnetic induction jumps AS'; and AS from the inter- 
section criterium of the full low- and high-temperature 
expressions (jSJ and {T)) without further approximation. 
To accomplish this, we use the elastic constant cn given 
by Brandt in Ref. [2ll |. The intersection criterium 
of the low and high temperature expansion of the 
partition function is then at least in the case of BSCCO 
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FIG. 2: Entropy jump per double layer per vortex ASi (first 
row) and jump of magnetic induction field AB (second row) 
at the melting transition. The experimental values for YBCO 
are from Ref. 28 for AS',, Ref Jp for AB by SQUID ex- 
periments (SQUID), and Ref. [S^ by torque measurements 
Torque). The experimental values for BSCCO are from Ref. 
by SQUID measurements (SQUID) of the magnetic field 



and Ref. |23by microscopic Hall sensors (Hall). 

a complicated integral equation via the dependence of 
C44 on the Lindemann parameter cl- One can solve this 
integral equation by numerical methods. The results are 
shown in Fig. 1 and Fig. 2. 



IV. DISCUSSION 

Our approximate analytic results H16|) and H19() for 
YBCO and BSCCO turn out to give practically the same 
curves. For comparison, we show in both figures the ex- 
perimental curves for YBCO of Ref. ,28^, .2£, ^ and for 
BSCCO of Ref. IH [H The good agreement in Fig. 1 
with the theoretical curves based on Eq. I|13|) shows that 
the Lindemann number is independent of the magnetic 
field for YBCO for the entire temperature range. For 
BSCCO the agreement is good for smaller i3-fields, where 
the second term in H18|l indroduces some dependence of 
the Lindemann on B, the largest near B = Q. There is 
some disagreement in Fig. 1 at larger B and in Fig. 2 at 
small B. This is not surprising since our vortex lattice 
model cannot be a good approximation in these regimes. 
At high B, the discrepancy comes mainly from Joseph- 
son decoupling of the layers 32], most pronounced for 
the strongly anisotropic BSCCO su per conductor, which 
leads also to large pinning effects [33. We think that 
this is also the reason for the difference in the curves of 
Kadowaki et al. in Ref. |^ and of Zeldov et al. in Ref. 
27 shown in Fig. 2. Pinning has the largest influence on 
the form of the melting curve at high B 34] , resulting in 
a decrease of AS; and AB in the limit of low temper- 
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atures shown by the curves of Zeldov et ai. Near B = 0, 
our model does not include the increase of the entropy 
coming from the thermal creation of vortices, in addition 
to the ones caused by the external magnetic field which 
form the lattice [s^. Moreover, in YBCO order parame- 
ter fluctuations become important |3fi| | which are ignored 
here. 

Summarizing, we have obtained the melting curve, the 
entropy, and the magnetic jump from a simple lattice 
defect model, and derived the Lindemann rule, including 
the size of the Lindemann number, and corrections to 
it. The determination of jump quantities over the phase 



transition cannot be obtained by the simple Lindemann 
rule. Our curves agree well with the experimental curves 
for YBCO and BSCCO except at zero and large magnetic 
fields. The simplicity of the model has allowed us to 
derive all results via analytic approximations. Our defect 
model is the simplest extension of the linear elasticity 
theory of vortex displacements. We have merely added 
integer- values defect gauge fields which introduce into the 
elasticity theory the rich physics of other phases of the 
vortex lattice caused by defect fluctuations, in particular 
the liquid phase and the associated melting transition. 
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